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ABSTRACT 

We analyze the relation between the mass of the central supermassive black hole {Mbh) and the 
number of globular clusters (Ncc) in elliptical galaxies and bulges as a ramification of the black hole 
fundamental plane, the theoretically predicted and observed multi-variable correlation between Mbh 
and bulge binding energy. Although the tightness of the Mbh-Nqc correlation suggests an unlikely 
causal link between supermassive black holes and globular clusters, such a correspondence can exhibit 
small scatter even if the physical relationship is indirect. We show that the relatively small scatter of 
the Mbh-Noc relation owes to the mutual residual correlation of Mbh and Nqc with stellar mass 
when the velocity dispersion is held fixed. Thus, present observations lend evidence for feedback- 
regulated models in which the bulge binding energy is most important; they do not necessarily imply 
any 'special' connection between globular clusters and Mbh- This raises the question of why Nqc 
traces the formation of ellipticals and bulges sufficiently well to be correlated with binding energy. 
Subject headings: black hole physics — galaxies: evolution — galaxies: formation — galaxies: elliptical 
and lenticular, cD 



1. INTRODUCTION 

There are now well-established correlations be- 
tween the mass of supermassive black holes (SMBHs) 
and properties of their host galaxies, such as bulge 
luminosity, mass, light concentration, and velocity dis- 
persion (IKo rmendv & Richstonc 1995; Magorrian ct al. 
19981: lE^rarese fc MerrittI 120001: IGebhardt et al.l 
2000t iTremaine et al.l l200a iMarconi fc Huntl 120031: 
Giiltekin et al.l |2009() . This suggests that the physical 



mechanism driving growth of the SMBH also plays a key 
role in forming the bulge (for spiral galaxies) or galax 



(for ellipticals) . Anal ytical estimates (Silk & Reei 
[Burkcrt fc Sil^ 120011 : Hopkins fc Hernauist 200^ 
well as numerical simulations (Di Mattco ct al 
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as 



Sprine^el et"all l2005t iCox et al.l ,2006,: .Robertson et al.l 
20061: ICroton et al.l [2006: Joha nsson et al.l I2009D with 
simple prescriptions for SMBH accretion have demon- 
strated the plausibility of this inference by matching the 
expected slopes of these correlations. 

Regardless of the detailed feedback prescription, these 
models predict that SMBHs grow until reaching some 
critical mass, where the energy and/or momentum re- 
leased by feedback expels material from the nucleus. As 
such, they robustly predict that the "true" correlation 
should be between SMBH mass and a quantity such as 
the binding e nergy or poten t ial wel l depth of material 
in the bulge. iHopkins et al.l ()2007al ) show that the ob- 
served correlations with different variables, and impor- 
tantly their scatter and systematic deviations from the 
relations, can be understood as the projections of a sin- 
gle fundamental dependence. This relation is approxi- 
mated closely by a multi- varia ble correlation, a bla c k hole 
fundamental plane (BHFP). lAller fc Richstond (|2007f ) 
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confirmed this in a sample of ellipticals and spiral 
bulges using dynamical models of bulge potentials, and 
iFteoli fc Mancini (2009) did so with simple proxies such 
as Mbh oc i?b ~ M^a^. 

Additional correlations have been found be- 
tween SMBH mass and dark matter halo mass, 
as well as the number Nqc of globular clusters 
(GC) in th e host galaxy fcp itler fc Forbes 200j; 
iBurkert fc Trcmai^ ^ I2OIOI: iHarris fc Harris! I2O1O0 . In 
particular, .Burkert fc Tremaind (|2010L hereafter BTIO) 
argued that Nqc is a better predictor of Mbh than 
the velocity dispersion ct, citing a smaller intrinsic 
scatter and a residual correlation between Nqc and 
Mbh in elliptical galaxies even after accounting for the 
median Mbh — c correlation, suggesting a fundamental 
link between the accretion of gas by the SMBH and 
the formation of a g alaxy's globular cluster system. 
IHarris fc Harris! ()2010L hereafter HHIO) extended the 
sample by making reasonable estimates of Nqc from 
the literature in galaxies with Mbh measurements. 

In this letter, we illustrate that the above link can be 
understood as a consequence of the BHFP relation com- 
bined with a residual correlation between Nqc and the 
bulge's stellar mass at fixed ct. Rather than suggest- 
ing a single "best" correlation between Mbh and a single 
galaxy parameter, the BHFP implies that the best pre- 
dictor of SMBH mass is some combination thereof. For 
example, Mbh has a positive correlation with the bulge's 
stellar mass even at fixed ct. Although the number of 
globular clusters in a particular galaxy, like Mbh, is a 
complex function of the galaxy's formation history, there 
exists a similar positive residual correlation between Nqc 
and M,, so that the resulting Nqc-Mbh residuals (fix- 
ing ct) will be positively correlated. 

In j |2] we descr i be a sample of 32 elliptical galaxies 
from Pen g et aD ( 20081) with auxiliary data compiled 
in Hopkin s et al.l (|2008bl and subsequent papers). In 
121 we fit separately the relations M^~a and Ncc"'^ in 
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these galaxies to establish the residual correlation be- 
tween Nqc and . Then we combine this residual slope 
with knowledge of the Mbh-M^ correlation at fixed a 
from the BHFP, and calculate the residual correlation 
and scatter expected between Nqc and Mbh- We sum- 
marize and conclude in Jj4l 

2. THE DATA 

To determine the dependence of Nqc on at fixed 
g, we cross-match objects compiled in iHopkins et alJ 
(|2008b( ) and subsequen t works w i th th e ACS Virgo 
Cluster Su r vey ( VCC, ICote et alJ [2001 . from which 
iPeng et all ()2008[ ) determined globular cluster counts 
{Ngc) and uncertainties. Following BTIO we obtained 
Ngc for several additional galaxies from iSpitler et alJ 

lioos). 

We obta i ned st ellar masses and un certainties from 
IPeng et all (|2008f l and IHopkins et all (2009a ,h). who 
compi led photome t ric d at a from several authors 
(e.g. iBender et alJ [19881 iRothberg fc JosephI l200l 
iLauer et al.l 120071: iKormendv et al.ll2009l and references 
therein). We use velocity dispersions as compiled by 
IHopkins et al.l (2009a,b). The latter quantity is the one 
best-determined for nearby massive galaxies, so we as- 
sume a log-uniform uncertainty in a of 0.02 dex, con- 
sistent with literature values. This approach yields 33 
galaxies for which we will determine the residual correla- 
tion between Nqc and . We discard the known recent 
merger remnant NGC1316 because its nuclear velocity 
dispersion is unrela xed and the glob ular cluster system 
is actively evolving (jSchweizeii ll 980f ) . The galaxy prop- 
erties used to analyze this 32-galaxy sample are provided 
in Table [TJ 

In addition, we will utilize the Mbh, Nqc, and a 
data directly from Table 1 of BTIO (compiled mostly 
from iGultekin et al.l ()2009f )). and the Mbh, Nqc data 
from Table 1 of HHIO. For the latter, we use the 
mean recorded a values froin the Hyperleda database 
(|Paturel et al.ll2003t IMcEhovl [19951) . These two sources 
yield 21 galaxies that serve as a combined comparison 
sample (Table[2]) for our derived Mbh^Nqc residual cor- 
relation. We note that BTIO and HHIO use slightly dif- 
ferent, but statistically consistent, values for Nqc where 
the samples overlap. The differing values of Mbh may 
make a larger difference in cases where multiple measure- 
ments exist; here, we follow BTIO and give half weight 
to each in our fits. Furthermore, different studies pro- 
vide different values for the velocity dispersion of a given 
galaxy; for example, Hyperleda returns a values ^^10-20 
km/s smaller than the ones from BTIO, and the papers 
by Hopkins et al. provide values that differ by ^ ±10-20 
km/s. We computed the observed residual Nqc^Mbh 
correlation using these alternate sources of a, and find 
that the small changes this introduces leave our conclu- 
sions completely unchanged. 

3. CORRELATIONS 

Both the BHFP and Mbh-Nqc relation appear to 
perform better than the Mbh~(^ relation because they 
explain its residuals and hence have a smaller intrinsic 
dispersion. Rigorously, this can be restated as follows: 
at fixed it, the residuals in A'Ibh correlate tightly with 
the residuals in Nqc and in for the BHFP. 

The question then arises: does Nqc directly explain 



the residuals in Mbh, or could the latter be attributed 
to other variables already proposed? Specifically, we ex- 
amine whether the observed Nqc-Mbh relation is pre- 
dicted as an indirect consequence of the BHFP relation. 
This is motivated by the BHFP prediction that at fixed a, 
good tracers of the bulge binding energy correlate tightly 
with the residuals in Mbh- (Interestingly, a similar re- 
lation was shown for the binding energies of individual 
Milkv Wav globulars bv iMcLaughliiil pOOOf ).) An ap- 
proximation to the bulge binding energy in ellipticals is 
a quantity like M^.a'^, so that at fixed a, galaxies with 
larger Af* will have a greater binding energy. Thus, if 
Nqc adequately traces at fixed a, as might be nat- 
ural given the observed A'^G'c~galaxy correlations, then 
the tightness of the Nqc-Mbh relation is expected. 

In this work, we focus on a particular projection of 
the BHFP that uses the bulge binding energy as the 
driving parameter. However, we note that the gen- 
eral BHFP, and also the corresponding sa l ient relation 
for G Cs (e.g. [Harris fc van den BCTghl[l98ll: [McLaughlinI 
[19991 and subsequent works) depends on galaxy forma- 
tion h istory in a more complicated way ([Hopkins et al.l 
l2009c[ ). Thus while the bulge binding energy serves ade- 
quately for our purposes, a more detailed accounting of, 
for example, the total baryon mass may lead to an even 
tighter expected correspondence. 

As follows, we calculate the expected residuals in Ncc~ 
Mbh assuming that this relation is a consequence of 
the BHFP and no other physics. As an expression 
of the BHFP correlation, we use the relation between 
the mass of the SM BH and bulge binding energy from 
IHopkins et all (|2007ai) . 

\ogMBH = r^ + p\og{M,a'^), 

where r/ = 8.23 ± 0.06, and P = 0.71 ± 0.06. We wih 
denote this quantity as predicted from the other ob- 
served variables as log(MBj:/|BHFP). Then we subtract 
from this the logarithm of Mbh as predicted solely from 
Mbh-o', denoted by \og{MBHW) , to obtain the BHFP- 
predicted residual correlation between Mbh and M^. We 
will signify this difference in logarithmic quantities as 
AMbh (AM,, ANgc)-- 

AMbh = \og{MBH\BHFP) - \og{MBH\<j) - 
= /31og(M,a2) - /31og((Af,|a)a2) 

= /3(log Af, - log(Af,|a)) = (3AM,. (1) 

This is just the statement that at fixed a, Mbh oc M^ . 
If there exists a relation between ANqc and AM,: 

ANgc ^1 AM,, (2) 

this will therefore result in a correlation between ANqc 
and AMbh- 

ANgc = ^AMbh - aAMnn- (3) 

To test this, we first calculate 7 from the existing data 
to establish a correlation between Ncc and M, at fixed 
velocity dispersion a (Figure [1]), and then combine it 
with the BHFP to create a prediction for the Ngc~Mbh 
residuals. Then, in Figure [H we compare this prediction 
to the observed residual correlation between Ncc and 
Mbh from BTIO and HHIO. 
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Figure 1. Top; Observed correlations of Nqc a-nd Af» witli velocity dispersion cr for the 32-galaxy sample described in JS] Lower Left: 
Correlation between Nqc and M* at fixed cr: AA''g'(7 (AM«) is the difference between the observed logarithm of Nqc (Af*) and the 
expected value of the logarithm of Nqc {Mt) given the linear relation in the top left (top right) panel. The lines are our regression curves 
fitted to the data. The dashed black curves are fitted using a technique accounting for intrinsic sca tter l|Tremaine et al. 2002), while the 
dashed red curves use an alternative maximum-likelihood technique from lAkritas fc Bershadvl (|1996l) . There is a clear positive correlation 
between ANqq and AM*. If the BHFP is the true underlying relation, then this observed residual correlation will lead to a correlation 
between Nqc and Mgjj, even at fixed a and with no other physics linking Mgu to the globular cluster systems. Lower Right: The 
observed correlation between Nqq and A/,. 



Specifically, we begin by determining the observed 
best-fit linear correlations between logA^cc and logcr, 
and logAf, and logcr. The resulting fits are shown in 
the top two panels of Figure [TJ We undertake all fits 
using two methods: the x^-minimization methods of 
iTremaine et al.l (|2002D (hereafter, T02), and the bivari- 
ate correl ated errors and intrinsic sc atter (BCES) esti- 
mators of lAkritas fc Bershadvl (|1996[ ). For the former, 
we account for intrinsic scatter in the Y axis by adding 
a uniform scatter in quadrature with the measurement 
errors such that the reduced-x^ value of the fit is 1. 
For the latter, we choose the regression line that bisects 
the BCES(Y|X) and BCES(X|Y) curves. Ahering these 
choices leads to small changes in the residual values, but 
does not change the residual slope in a statistically sig- 
nificant way. One-sigma uncertainties in the fitted slopes 
are calcu lated using paired n onparametric bootstrap sim- 
ulations (jBabu fc Raolll993[ ). 

For each object, we then use the fitted relation to com- 
pute the expected value of log A^gc (logM*) given its 
observed value of cr, and subtract it from the observed 
value of log A^GC (logM*) to obtain AA^gc (AM*). In 



the lower left panel of Figure[I] as expected we see a clear 
positive correlation between A A^gc and AM* , indicating 
that at a fixed cr, elliptical galaxies with more globular 
clusters also have a larger total stellar mass. We note 
that the values of ANqc and AMbh calculated using 
BCES or x^-niinimization on the direct correlations are 
the same to within 0.1 dex. Subsequent estimates of the 
residual slope are unaffected by this choice. 

However, the two fitting methods obtain somewhat dif- 
ferent estimates for the slope of the resulting residual cor- 
relation, 7, which we will use to estimate the expected 
Nqc^Mbh residual slope a = j/p. As plotted in Fig- 
ure [TJ we find 

7j,02 = 0.79 ± 0.25 
IBCES - 1.11 ±0.21, 
and correspondingly, 

aro2 = 1.11 ±0.36 
ascES = 1.56 ±0.32. 

In Figure [21 we compare this predicted slope with the 
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Figure 2. Residual correlation between Nqc and Mbh at fixed velocity dispersion a. The points are data from BTIO and HHIO; we 
follow BTIO where multiple Mgjj measurements exist by assigning each measurement half weight in any fits. Red points correspond to 
the BTIO sample; their positions here and in BTlO's Figure 3 differ slightly because we take Nqc residuals with respect to a instead of 
bulge luminosity, but the residuals here have roughly the same slope and span a similar range of ANqc ^ compared with BTIO. The 
eight points added by the elliptical sample of HHf reinforce this residual correlation and expand its range. The dark gray shaded region 
is the predicted residual correlation (with associated 1-sigma uncertainty in lighter shade) assuming that Mgjj is determined only by the 
BHFP, combined with the observed correlation between Nqc and M* at fixed tr determined by fitting the data in Figure [T] We see that 
the observed residual slope is in good agreement with this expected slope. Thus the apparent additional predictive power of Nqc for ^BH 
can be entirely accounted for by the predicted correlation of Mbh and the observed correlation of Nqc with the bulge binging energy. 



one observed in the comparison sample of BTIO and 
HHIO. Again, we compute the residuals against a in 
both Ngc and Mbh', note that for consistency, this is 
slightly different than the quantities plotted in BTlO's 
Figure 3 where the Nqc residual was computed against 
the bulge luminosity, not velocity dispersion. We plot 
a region in light gray to highlight the extremes of the 
predicted slopes, corresponding to the range bounded by 
the 1-sigma uncertainties in the slope given by our two 
regression methods. In darker gray we simply plot the 
range bounded by our two slope estimates. 

As in Figure [T] we fit the data directly and find that 
the observed residual slope between Nqc and M bh is 

arm = 0.78 ± 0.24 
ctBCES = 1-33 ± 0.34, 

in good agreement with the BHFP predictions above. 

The detection of this residual correlation by BTIO 
quantitatively demonstrates that Nqc is a better pre- 
dictor of Mbh than is a. Such a comparison can be 
alternatively phrased as a reduction in the intrinsic scat- 
ter of the correlation. In BTIO, the intrinsic scatter of 
Mbh~<^ was found to be e ^ 0.3 dex, while the intrinsic 
scatter in Nqc-Mbh is e ~ 0.2 dex. The magnitude of 



this dispersion can be predicted by combining the BHFP 
relation with the obse rved correlations with Nqc- From 
ms "eraII(|2007aD . we see that the Mbh~Ei, correla- 
tion has an intrinsic scatter ~ 0.2-0.25 dex, and from the 
present data, we find the scatter of Et-Ncc is 0.22 ±0.04 
dex. By propagating these as measurement uncertainties 
to the Mbh-Ngc relation, we predict that the measured 
Nqc-Mbh intrinsic scatter should be e = 0.23 ± 0.03 
dex, consistent with the measurement by BTIO. It is 
also consistent with the combined dataset of BTIO and 
HHIO, for which we find e = 0.21 ± 0.04 dex. 

4. CONCLUSIONS 

We have shown that the number of globular clusters 
in elliptical galaxies exhibits a residual dependence on 
at fixed tr, implying that the bulge binding energy 
(~ MH.tT^) is a better indicator of Nqc than a or 
alone. The same wa s shown to be true for Mbh by 
iHopkins et al.l (|2007al ). as these parameters constitute 
a formulation of the BHFP. Thus the apparent power 
of Mbh-Ngc versus Mbh-o' owes to the fact that Nqc 
and Mbh are both tracers of the same fundamental prop- 
erty such as the bulge binding energy. 

This resolves several puzzling aspects of the previous 
interpretation of the data. As BTIO themselves point 
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out, there cannot be a direct causal correlation between 
Ngc and Mbh, since most of the GC mass is at very 
large radii and has never had any interaction with the 
galaxy nucleus. Moreover, while most GCs likely formed 
at very high redshift, the final mass of the SMBH is 
sensitive to its gro wth via gas accretion at z S, 2 (e.g . 
iHopkins fc Hernaui st 2006; Hopkins et al.ll2007cri2008a[ ). 
However, this naturally predicts that Nqc should serve 
reasonably well as a mass tracer, so that the depen- 
dence of Mbh on and formation time leads to a 
surprisingly tight but expected Nqc-M bh correlation. 
The sa me arguments explain the result in IHopkins et all 
(|2009c[ ). who show that the observed Mbh is sensitive 
to the entire galaxy baryonic mass ~ i.e. perhaps the 
mass traced by Nqc is the same as the mass that ac- 
tually sets the escape velocity and potential well depth 
at i? = 0, rather than just the stellar mass enclosed in 
a small radius around the BH, which can vary widely in 
systems of similar Mbh- Such a relation between Nqc 
and global g alaxy mass or luminosity has been demon- 
strat ed (e.g. [Harris &: van den Bergh|[l981t iMcLaughlinI 
Il999f) . and this trait supports the idea that Nqc and 
Mbh are connected indirectly by a more fundamental 
galaxy property. 

This also naturally explains why HHIO find that the 
relation breaks down for SO galaxies. These galaxies are 
structurally different than ellip ticals and may have dif- 
ferent formation histories (Lars on et al.i ri980'). so Nqc 
and the total stellar mass may not faithfully trace the 
bulge binding energy. Sin ce SO's are not partic ularly dis- 
crepant in Mbh-c (e.g. iGiiltekin et al.l l2009l and pre- 
vious works), this suggests that the TVcc-bulge relation 
is the connection that weakens for these systems. HHIO 
also find no statistically significant correlation in spirals: 
although three out of the four spirals from HHIO lie on 
the Ngc~Mbh relation, there simply isn't yet enough 
data to know for sure if this relation persists for spi- 
ral bulges. However, the underlying BHFP relation ties 
Mbh to the binding energy and explains its residual 
correlations with bulge parameters, even for these disky 
galaxies where Nqc possibly deviates. This alone sug- 
gests that the BHFP, not Mbh-Nqc or Mbh-<^, is the 
' more fundamen t al' corr elation. 

IHopkins et al.l ()2007b[) showed that the existence of a 
black hole fundamental plane is a robust prediction of nu- 
merical simulations of gas-rich mergers that include the 
effects of gas dissipation, cooling, star formation, and 
black hole accretion and feedback. The present work 
shows that this local and widely expected correlation be- 
tween supermassive black hole mass and bulge binding 
energy in feedback-regulated scenarios, combined with 
a similar correlation for NgCi can account for the ob- 
served Ngc^Mbh relation and its scatter. The inter- 
esting question raised by such a correlation is not why 
Ngc correlates tightly with Mbh, since this is indirect, 
but why Ngc correlates tightly with galaxy binding en- 
ergy/potential well d epth. Some such correlation is ex- 
pected and observed (|McLaughlinlll999l: IBlakesled 11991 
IPeng et "al. 2008): an example is that systems at fixed ve- 
locity dispersion with higher stellar mass have accreted or 
formed more stars, likely including globular clusters. But 
that the iVGc-bulge relation should be so tight, and in- 
clude both metal-rich and metal-poor populations, may 
support the inferences by BTIO and HHIO (and refer- 



ences therein) that the formation of globular cluster sys- 
tems and growth of supermassive black holes in elliptical 
galaxies are driven by a common galaxy property. 
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comments. 

REFERENCES 

Akritas, M. G., & Bershady, M. A. 1996, ApJ, 470, 706 

AUer, M. C, & Richstone, D. O. 2007, ApJ, 665, 120 

Babu, G. J., & Rao, C. R. 1993, Handbook of Statistics, Vol. 9 

(Amsterdam: Elsevier Science), 627-659 
Babu, G. J., & Singh, K. 1983, Annals of Statistics, 11, 999 
Bender, R., Doebereiner, S., & Moellenhoff, C. 1988, A&AS, 74, 

385 

Blakeslee, J. P. 1999, AJ, 118, 1506 
Burkert, A., & Silk, J. 2001, ApJ, 554, L151 
Burkert, A., & Tremaine, S. 2010, ApJ, 720, 516 
Cote, P., et al. 2004, ApJS, 153, 223 

Cox, T. J., Dutta, S. N., Di Matteo, T., Hernquist, L., Hopkins, 

P. F., Robertson, B., & Springel, V. 2006, ApJ, 650, 791 
Croton, D. J., et al. 2006, MNRAS, 365, 11 

Di Matteo, T., Springel, V., & Hernquist, L. 2005, Nature, 433, 
604 

Feoli, A., & Mancini, L. 2009, ApJ, 703, 1502 
Ferrarese, L., & Merritt, D. 2000, ApJ, 539, L9 
Gebhardt, K., et al. 2000, ApJ, 539, L13 

Giiltekin, K., et al. 2009, ApJ, 698, 198 

Harris, G. L. H., & Harris, W. E. 2010, 'arXiv:1008.4748j 
Harris, W. E., & van den Bergh, S. 1981, AJ, 86, 1627 
Hopkins, P. F., Cox, T. J., Dutta, S. N., Hernquist, L., 

Kormendy, J., & Lauer, T. R. 2009a, ApJS, 181, 135 
Hopkins, P. F., Cox, T. J., Keres, D., & Hernquist, L. 2008a, 

ApJS, 175, 390 
Hopkins, P. F., & Hernquist, L. 2006, ApJS, 166, 1 
Hopkins, P. F., Hernquist, L., Cox, T. J., Dutta, S. N., & 

Rothberg, B. 2008b, ApJ, 679, 156 
Hopkins, P. F., Hernquist, L., Cox, T. J., Robertson, B., & 

Krause, E. 2007a, ApJ, 669, 45 
— . 2007b, ApJ, 669, 67 

Hopkins, P. F., Lauer, T. R., Cox, T. J., Hernquist, L., & 

Kormendy, J. 2009b, ApJS, 181, 486 
Hopkins, P. F., Murray, N., & Thompson, T. A. 2009c, MNRAS, 

398, 303 

Hopkins, P. F., Richards, G. T., & Hernquist, L. 2007c, ApJ, 654, 
731 

Johansson, P. H., Naab, T., & Burkert, A. 2009, ApJ, 690, 802 
Kormendy, J., Fisher, D. B., Cornell, M. E., & Bender, R. 2009, 

ApJS, 182, 216 
Kormendy, J., & Richstone, D. 1995, ARA&A, 33, 581 
Kundu, A., & Whitmore, B. C. 2001, AJ, 121, 2950 
Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, ApJ, 237, 

692 

Lauer, T. R., et al. 2007, ApJ, 664, 226 
Magorrian, J., et al. 1998, AJ, 115, 2285 
Marconi, A., & Hunt, L. K. 2003, ApJ, 589, L21 
McElroy, D. B. 1995, ApJS, 100, 105 
McLaughlin, D. E. 1999, AJ, 117, 2398 
— . 2000, ApJ, 539, 618 

Paturel, G., Petit, C., Prugniel, P., Theureau, G., Rousseau, J., 
Brouty, M., Dubois, P., & Cambresy, L. 2003, A&A, 412, 45 
Peng, E. W., et al. 2008, ApJ, 681, 197 

Robertson, B., Hernquist, L., Cox, T. J., Di Matteo, T., Hopkins, 
P. F., Martini, P., & Springel, V. 2006, ApJ, 641, 90 

Rothberg, B., & Joseph, R. D. 2004, AJ, 128, 2098 

Schweizer, F. 1980, ApJ, 237, 303 

Silk, J., & Rees, M. J. 1998, A&A, 331, LI 

Spitler, L. R., & Forbes, D. A. 2009, MNRAS, 392, LI 

Spitler, L. R., Forbes, D. A., Strader, J., Brodie, J. P., & 
Gallagher, J. S. 2008, MNRAS, 385, 361 

Springel, V., Di Matteo, T., & Hernquist, L. 2005, MNRAS, 361, 
776 

Tremaine, S., et al. 2002, ApJ, 574, 740 



Table 1 

Galaxy properties for Figure [Tl 



Galaxy 


Ngc 


Mt 


(7 








(km/s) 


NGGU821 


320 ± 45 


229 ± 57 


209 ± 10 


NGGld99 


5800 ± 700 


363 ± 91 


359 ± 18 


NGC3377 


266 ± 66*" 


26.3 ± 6.6 


141 ± 7 


NGC3379 


270 ± 68=^ 


107 ± 27 


221 ± 11 


NGC4318 


18 ± 6.1 


5.0 ± 1.7 


101 ± 5 


NGC4365 


3246 ± 598 


226 ± 52 


269 ± 13 


NGC4374 


4301 ± 1201 


236 ± 61 


287 ± 14 


NGC4382 


1110 ± 181 


186 ± 44 


196 ± 10 


NGC4387 


69.5 ± 9.8 


13.7 ± 3 


84 ± 4.2 


NGC4406 


2660 ± 129 


289 ± 60 


250 ± 12 


NGC4434 


141 ± 34 


21.4 ± 4 


118 ± 6 


NGC4458 


72 ± 12 


8.7 ± 2 


85 ± 4.3 


NGC4459 


218 ± 28 


77.9 ± 14 


168 ± 8 


NGC4464 


25.3 ± 9.2 


7.1 ± 1.4 


120 ± 6 


NGC4467 


-6 ± 13 


1.8 ± 0.6 


67 ± 3.4 


NGC4472 


7813 ± 830 


531 ± 110 


287 ± 14 


NGC4473 


76 ± 97 


53.5 ± 12 


178 ± 9 


NGC4476 


20.1 ± 7.3 


3.7 ± 2 


41 ± 2.1 


NGC4478 


58 ± 11 


22 ± 4 


149 ± 7 


NGC4486 


14660 ± 891 


302 ± 79 


360 ± 18 


NGC4489 


31 ± 9 


6.98 ± 1.7 


49 ± 2.4 


NGC4515 


81 ± 10 


7.7 ± 1.5 


90 ± 4.5 


NGC4551 


47 ± 11 


11.6 ± 2.4 


100 ± 5 


NGC4552 


984 ± 198 


95 ± 16.9 


261 ± 13 


NGC4564 


213 ± 31 


26 ± 6 


153 ± 8 


NGC4621 


803 ± 355 


83.9 ± 19 


237 ± 12 


NGC4649 


4745 ± 1099 


339 ± 50 


341 ± 17 


NGC4660 


205 ± 28 


23.8 ± 4 


191 ± 9 


VCC1199 


-9 ± 14 


0.58 ± 0.16 


69 ± 3.5 


VCC1440 


26.7 ± 6.8 


1.2 ± 0.44 


59 ± 3 


VCC1627 


3.6 ± 3.7 


1.0 ± 0.32 


47 ± 2.4 


VCC1871 


10.4 ± 5 


2.3 ± 0.58 


51 ± 2.6 



Note. — Properties of elliptical galaxies used to 
determine the correlation of Nqc with M, at fixed 
velocity dispersion a (Figur e [T]|. Values of Nqc 
and Mt are compiled from IPeng et al.l 1)20081 ) un- 
less otherwise noted. Values of a are compiled from 
IHopkins et al.l ll2009a| [hD and assumed to have a log- 
uniform uncertainty of 0.02 dex. 

"Spitlcr et alT (lOOS") 
fa .Kundu fc Whitmorj l[200ll 'l 



Globular Cluster - Black Hole Relation 



Table 2 

Galaxy properties for Figure [2] 



Galaxy 


Mbh 


-f 1-sigma 


-1-sigma 


Nnn 


(J 


ANnn 






(Mo) 




(Mo) 




(km/s) 








Tremaine (2010) 












NGC0821 


4.2 X lO'^ 


2.8 X 10"^ 


8 X 10° 


320 ± 45 


209 lb 10 


-0.39 


-0.79 


NGC1316 


1.5 X 10* 


8 X 10^ 


8 X lO'^ 


1173 ± 240 


226 ± 9 


0.05 


-0.37 


NGC1399 


1.3 X 109 


5 X 10^ 


7 X 10* 


5800 ± 700 


337 zb 16 


0.11 


-0.13 




5.1 X 10* 


/ X iU 


7 X lO'^ 






0.11 


-0.54 


NGC3377 


1.1 X 10* 


1 1 V 1 

1 . ± X ±u 


1 X lO'^ 


266 ± 66 


145 ± 7 


0.11 


0.27 


IN vjv_jOO I y 


1.2 X 10* 


O X iU 








-0 44 


-U.Ol 




4 X 10* 


1 X ±u 


1 X 10* 






-0.44 


0.22 


NGC4374 


1.5 X 10^ 


T 1 V 1 

1 . i X iU 


6 X 10* 


4301 ± 1201 


296 ± 14 


0.19 


0.16 


NGC4459 


7.4 X lO'^ 


T /I V 1 n7 

1.4 X iU 


1.4 X lO'^ 


218 ± 28 


167 ± 8 


-0.20 


-0.15 


NGC4472 


1.8 X 10^ 


X iU 


6 X 10* 


7813 ± 830 


310 ± 10 


0.37 


0.16 


NGC4486 


6.4 X 10^ 


5 X 10* 


5 X 10* 


14660 ± 891 


375 ± 18 


0.34 


0.38 


NGC4564 


6.9 X lO'^ 


4 X 10^ 


1 X lO'^ 


213 ± 31 


162 ± 8 


-0.17 


-0.13 


NGC4594 


5.5 X 10* 


5 X 10^ 


5 X lO'^ 


1900 ± 189 


240 ± 12 


0.16 


0.09 


NGC4649 


4.5 X 10^ 


1 X 10^ 


1 X 10^ 


4745 ± 1099 


385 ± 19 


-0.19 


0.18 


NGC5128 


3 X 10* 


4 X 10^ 


2 X lO'^ 


1550 ± 390 


150 ± 7 


0.82 


0.65 




7 X lO'^ 


1.3 X 10^ 


3.8 X lO'^ 






0.82 


0.01 


Harris & Harris (2010) 














NGC2778 


1.6 X lO'^ 


9 X 10*5 


2 X 10^ 


50 ± 30 


162 ± 8 


-0.80 


-0.76 


NGC4261 


5.5 X 10* 


1.1 X 10* 


1.2 X 10* 


530 ± 100 


309 ± 14 


-0.79 


-0.35 


NGC4473 


1.3 X 10* 


5 X 10^ 


9.4 X lO'^ 


376 ± 97 


180 ± 8 


-0.08 


-0.03 


NGC4552 


4.8 X 10* 


8 X 10^ 


8 X lO'^ 


1200 ± 250 


253 ± 12 


-0.12 


-0.06 


NGC4621 


4 X 10* 


6 X 10^ 


6 X lO'^ 


800 ± 355 


225 ± 11 


-0.11 


-0.06 


NGC4697 


2 X 10* 


2 X 10^ 


2 X lO'^ 


229 ± 50 


171 ± 8 


-0.22 


0.24 


NGC5813 


7 X 10* 


1.1 X 10* 


1.1 X 10* 


1650 ± 400 


237 ± 11 


0.12 


0.22 


NGC5846 


1.1 X 10^ 


2 X 10* 


2 X 10* 


4700 ± 1200 


239 ± 11 


0.56 


0.40 



Note. — Properties of elliptical galaxies used to calculate the correlation of Nqc with 
Mbh at fixed velocity dispersion cr (Figure[2ll, following BTIO. Data for Nqc a-^d Mbh are 
compiled from BTIO and HHIO. Values of cr are as presented in BTIO for those galaxies, and 
values of a for the HH IO ellipticals are taken as the mean value recorded in the HyperLeda 
database UPaturel et al..,200awMcElrovia995i ). 



